UPRF-03-002 



Subtraction terms at NNLO 



Stefan Weinzierl 1 



Dipartimento di Fisica, Universita di Parma, 
INFN Gruppo Collegato di Parma, 43100 Parma, Italy 



Abstract 

Perturbative calculations at next-to-next-to-leading order for multi-particle final 
states require a method to cancel infrared singularities. I discuss the subtraction 
method at NNLO. As a concrete example I consider the leading-colour contributions 
to e + e~ — > 2 jets. This is the simplest example which exhibits all essential features. 
For this example, explicit subtraction terms are given, which approximate the four- 
parton and three-parton final states in all double and single unresolved limits, such 
that the subtracted matrix elements can be integrated numerically. 
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1 Introduction 



Present and future collider experiments will provide a large sample of multi-particle final 
states. In order to extract information from this data, precise theoretical calculations 
are necessary. This implies to extend perturbative calculations from next-to-leading order 
(NLO) to next-to- next-to-leading order (NNLO). Due to a large variety of interesting jet 
observables it is desirable not to perform this calculation for a specific observable, but to set 
up a computer program, which yields predictions for any infra-red safe observable relevant 
to the process under consideration. Compared to already existing NNLO calculations for 
specific (inclusive) observables like Drell-Yan or inclusive Higgs production [l]-[5], this 
requires to work with fully differential cross-sections. This is the major complication of 
the project. The three major challenges to accomplish this task are: the calculation of 
two-loop amplitudes, a method for the cancellation of infrared divergences and stable and 
efficient Monte Carlo techniques. 

The last three years have witnessed a tremendous progress in the calculation of two- 
loop amplitudes [6] -[9]. In particular the two-loop amplitudes for e + e~ — > 3 jets have been 
calculated [10, 11]. 

Up to now, less is known for the cancellation of infrared divergences at NNLO. Loop 
amplitudes, calculated in dimensional regularization, have explicit poles in the dimensional 
regularization parameter e = 2—D/2, arising from infrared singularities. These poles cancel 
with similar poles arising from amplitudes with additional partons, when integrated over 
phase space regions where two (or more) partons become "close" to each other. However, 
the cancellation occurs only after the integration over the unresolved phase space has been 
performed and prevents thus a naive Monte Carlo approach for a fully exclusive calculation. 
It is therefore necessary to cancel first analytically all infrared divergences and to use 
Monte Carlo methods only after this step has been performed. Infrared divergences occur 
already at next-to-leading order. In this case, general methods to handle the problem are 
known. Examples are the phase-space slicing method [12, 13] and the subtraction method 
[14]-[18]. As already mentionend, similar methods at NNLO are not yet known. The 
major difficulty arises from so-called double unresolved configurations, where three partons 
become degenerate simultaneously. A phase-space slicing approach to double unresolved 
configurations has been used in the calculation of the photon + 1-jet rate in electron- 
positron annihilation [19]. The calculational complexity of this process is somewhere in 
between a NLO calculation and a full NNLO calculation, as the relevant two-loop amplitude 
vanishes. 

In this paper I discuss the general setup for the subtraction method at NNLO. To 
illustrate the method, I use as an example the leading-colour NNLO contributions to 
e + e~ — > 2 jets. This is the simplest example which exhibits all essential features. For 
this example I give all necessary subtraction terms. These subtraction terms approximate 
the leading-colour Born contribution of e + e~ — > qggq in all double and single unresolved 
limits, and the leading-colour one-loop contribution of e + e~ — > qgq in all single unresolved 
limits. The subtracted matrix elements can therefore be integrated numerically over the 
appropriate phase space. Valuable information for the construction of the subtraction 
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terms is provided by the known behaviour of the tree and one-loop amplitudes in singular 
limits. The double unresolved limits of tree amplitudes have been considered in [19]-[26]. 
Single unresolved limits of one- loop amplitudes have been considered in [27]- [33]. However, 
for the subtraction method this information is not yet sufficient. (It is sufficient for a phase 
space slicing approach.) Since the subtraction terms are subtracted over the complete phase 
space, they must have the correct singular behaviour not only in the double unresolved 
cases, but also in the single unresolved ones. This requires to extend the calculations of 
the singular behaviour of tree amplitudes and to include also subleading single unresolved 
singularities. 

The method presented here is general and not restricted to the example of e + e~ — > 
2 jets. Subtraction terms for other splittings (like g — > ggg) and other kinematical con- 
figuartions (e.g. with partons in the initial state) can be worked out along the same lines. 

The subtraction terms still need to be integrated analytically over the appropriate un- 
resolved phase space. Although some of the occuring integrals are quite involved, recently 
developped integration techniques [34, 35] seem to indicate that these integrals can be done 
analytically. 

The knowledge of the subtraction terms will also be useful for an improvement of parton 
showering algorithms in event generators [36]. 

This paper is organized as follows: In the following section I discuss the general setup 
of the subtraction method at NNLO. Sect. 3 gives the relevant amplitudes for the example 
e + e" — > 2 jets. Sect. 4 reviews the singular limits of amplitudes. In sect. 5 I present the 
NNLO subtraction terms and in sect. 6 I discuss the mechanism of cancellations between 
the various terms. Finally sect. 7 contains the conclusions and an outlook. 



2 General idea 

In this section I outline the general setup for the subtraction method at NNLO. The NNLO 
contribution to an observable O whose LO contribution depends on n partons is given by 

{(D) NNLO = (1) 

J O n+2 (pi,...,p n +2) da ( nl 2 + J O n+1 (p l ,...,p n+1 ) d(j^ ] +1 + J O n (pi,...,p„) da^\ 

The various contributions are symbolically shown in fig. 1. The observable O has to be 
infrared safe. This requires that whenever a n + 1 parton configuration pi,...,p n +i becomes 
kinematically degenerate with a n parton configuration p' Xr --,p' n we must have 

O n+1 ( Pl ,..., Pn+1 ) -> n (pi, -X). (2) 

In addition, we must have in the double unresolved case (e.g. when a n + 2 parton configu- 
ration pi,...,p n +2 becomes kinematically degenerate with a n parton configuration p'i,...,p' n ) 

O n+2 (pi,...,p n +2) -> O n (p[,...,p' n ). (3) 
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Figure 1: Modelling of jets at next-to-next-to-leading order. The cone labeled with y cut 
represents the experimental cuts. At next-to-next-to-leading order a jet is modelled either 
by one, two or three partons. 

The subscript n for O n and da^ indicates that this contribution is integrated over a n- 
parton phase space. The various contributions d<j$ are given by 

j_(0) ( AO) * AO) \ ji 

' a(0) \ j. 



da W _ (AP) * .(i) , A W *A0) " 

d(7^ = A.^ + + «^n^ *^h) J d(j) n , (4) 

where .4*? denotes an amplitude with n external partons and I loops. d<p n is the phase 
space measure for n partons. Taken separately, each of the individual contributions da^ 2 , 

d&n+i an d dan" 1 gives a divergent contribution. Only the sum of all contributions is finite. 
However, since one would like to perform the integration numerically, one must first cancel 
analytically the divergences, before numerical integration over the different phase spaces 
can be envisaged. This is not a new phenomena at NNLO, it already occurs at NLO. We 
can consider an observable, whose LO contribution depends onn+1 partons. The NLO 
contribution is given by 

(0>X° = J 'o n+2 (p 1 ,...,p n+2 )da% + J 'o n +i(pi,-,Pn+i)^2i- (5) 

To render the two contributions finite, one adds and subtracts a suitable term [15] 

J O n+2 da® 2 + J O n+2 da^ = 

J (0 n+2 da% - O n+1 o dagg) + J (a+i da^l + O n+1 o daft?) ■ (6) 

The approximation dacn+x has to fulfill the following requirements: First, da^+l must be 
a proper approximation of da^_ 2 such as to have the same pointwise singular behaviour 
in D dimensions as da^ 2 itself. Thus, dcg+l acts as a local counterterm for da^} 2 an d 



one can safely perform the limit e — > 0. Secondly, da n ^ must be analytically integrable in 
D dimensions over the one-parton subspace leading to soft and collinear divergences. The 
integral 

fdoSH (7) 

contains then explicit poles in 1/e 2 and l/e which exactly cancel the poles from the loop 
amplitude in da^l v Note that the observable O in the approximation term is evaluated 
with a n + 1 parton configuration. This allows the integration over the one-parton subspace 
to be performed analytically, but requires a mapping of four- momenta from the n+2 parton 
configuration to the n + 1 parton configuration. This mapping has to satisfy momentum 
conservation and also has to keep all partons massless. In addition it has to have the right 
behaviour in the singular limits. Several choices for such a mapping exist [15, 37]. The 
notation O n +i o da^+l is a reminder, that in general the approximation is a sum of terms 

O n+1 oda { ^l = ^O^daH (8) 

and the mapping used to relate the n + 2 parton configuration to a n + 1 parton configu- 
ration differs in general for each summand. 

In a similar way, I rewrite the NNLO contribution as 

{0) N n NL ° = j \o n+2 da n %-O n+1 odatil-O n oda^) 

+ J (a+i dali, + O n+l o da ( Xl - O n o da^) 

+ j (O n da™ + O n o da^ + O n o da^) . (9) 

Here dan' 2 ^ acts as a local counterterm for double unresolved configurations to da^+2 ~ 
da^+i in D dimensions. Furthermore, it is integrable over a two-parton subspace. dan' 1 ^ is 
a local counterterm to da^ 1 ^ +da n d +i in D dimensions, whenever one of the (n + 1) partons 
becomes unresolved. It is integrable over a one-parton subspace. This setup extends to 
higher orders and in general, dati'^ is an approximation to 

da {l) + da m + da (1 ' ,_1) + + da (,_1,1) - da (,,1) - da {l ' 2) - - da {l ' k ~ l) (10) 

e.g. a contribution with at most n + k partons and at most I loops in the case when k 
partons become unresolved. 

The task at NNLO is to give the explicit forms of da>n' 2 ^ (and of dai 0,1 * 1 and da n 1,V> ) and 
to integrate these approximations over a two-parton phase space (or a one-parton phase 
space in the case of dan' 1 ^ and dan' 1 ^). 
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To illustrate the cancellations for double unresolved contributions, I consider a simple 
toy model, where the double real emission contribution to a 1-jet observable is given by 

3 daf ] = 9 3 (x 1 ,x 2 ,x 3 )F 3 (x 1 ,x 2 ,x 3 ). (11) 

The singular behaviour of the functions F 3 is assumed to be 

1 

lim F 3 (x 1 ,x 2 ,x 3 ) = —F 2 (x 2 ,x 3 ), 

xi->0 X\ 

lim F 2 {x u x 2 ) = —F 1 (x 2 ). (12) 

xi-»0 X\ 

The function 9 3) defining the observable satisfies 

lim 9 3 (x 1 ,x 2 ,x 3 ) = 9 2 (x 2 ,x 3 ), 

lim 9 2 (x l ,x 2 ) = 9 1 (x 2 ). (13) 

Similar relations are assumed to hold for the singular limits in the other variables. The 
NLO subtraction term is given by 

2 o daf' l) = 

—F 2 (x 2 , x 3 )9 2 (x 2 , x 3 ) + — F 2 (x 1: x 3 )9 2 (x u x 3 ) + — F 2 (x x , x 2 )6 2 ( Xl , x 2 ). (14) 
Xi x 2 x 3 

Let us first consider the cancellation of single unresolved contributions to a 2-jet observable. 
I consider the x± — > limit. Here, 3 daf^ — 2 o daf' 1 ^ consists of two contributions: 
The first part 

lim ( F 3 (x 1 ,x 2 ,x 3 )9 3 (x 1 ,x 2 ,x 3 ) - — F 2 (x 2 ,x 3 )9 2 (x 2 ,x 3 ) J (15) 
is finite by construction. The second part is given by 

lim ( — -F 2 (xi,x 3 )6 2 (xi,x 3 ) - —F 2 (x 1 ,x 2 )6 2 (x 1 ,x 2 ) 

xi^O \ X 2 X 3 

1 -F 1 (x 3 )6 1 (x 3 ) - J—F 1 (x 2 )6 1 (x 2 ). (16) 



X ± X 2 

This part is proportional to 9i{x 2 ) or 9\{x 3 ) and therefore does not contribute to a 2-jet 
observable (9\ is zero for a two-jet observable). However, this argument no longer holds 
if we move to NNLO and consider double unresolved contributions to a 1-jet observable. 
Here, 9\ does not vanish and the r.h.s of eq. (16) gives a divergent contribution already 
in the single unresolved case x 1 — > 0, and x 2 ,x 3 7^ 0. At NNLO these singular pieces 
have to be taken properly into account and subtracted out. To summarize, the major 
complication for the subtraction method is not to find subtraction terms which match the 
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matrix element in all double unresolved limits, but to find subtraction terms, which match 
in all double unresolved limits and all single unresolved limits. For the toy example the 
NNLO subtraction term is given by 

O x o dap 2) = 

— F^e^xs) - —F 1 {x 2 )6 1 {x 2 ) - — iMxOflifci). (17) 

X X X 2 XxXz x 2 x 3 

Note that in this toy example the NNLO subtraction term can be obtained solely from the 
singular behaviour of the NLO subtraction term. In general however, the limits do not 
commute, e.g. taking the singular limit of the single unresolved approximation does not 
equal the double unresolved approximation of the matrix element. 



3 An example 

To make the method for the double unresolved contributions as transparent as possible, 
I consider as a simple example the leading iV c -contribution to the tree-level amplitude 
e + e~ — > qggq. This example is complicated enough to expose all essential features. 



3.1 Double real emission 

The singular behaviour of the kinematical structure of QCD amplitudes is entangled with 
colour factors. For a clear understanding of the factorization properties in soft and collinear 
limits it is desirable to disentangle the colour factors from the rest. This can be done either 
by introducing colour charge operators [15] or by decomposing the full amplitude into par- 
tial and primitive amplitudes. In this paper I follow the second approach, since this method 
is closer related to the way several one-loop amplitudes were calculated. Throughout this 
paper I use the following normalization for the colour matrices: 

Tr T a T b = ^5 ab . (18) 

Tree-level amplitudes are decomposed into colour factors and partial amplitudes. Each 
partial amplitude is gauge-invariant and has usually a simpler structure. For example, the 
colour decompostion of the the tree amplitudes for e + e~ — > q\g 2 g 3 q^ amplitude reads 



•AT (Ql,92, 93, Q4) = 

\T 2 T 3 ) u Af\ qu g 2 ,g 3 ,q 4 ) + (T 3 T 2 ) u A { °\ qi ,g 3 ,g 2 ,q 4 )} . (19) 



2 2 
e g c 



Here, c denotes a factor from the electro-magnetic coupling: 

c = -Q q + v e viV z (s), 
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The electron - positron pair can either annihilate into a photon or a Z-boson. The first 
term in the expression for cq corresponds to an intermediate photon, whereas the last term 
corresponds to a Z-boson. The left- and right-handed couplings of the Z-boson to fermions 
are given by 



vl 



h-Q sin 2 



w 



—Q sin 9w 



(21) 



sin 9w cos 9w cos 9w 

where Q and I3 are the charge and the third component of the weak isospin of the fermion. 
The short-hand notation for the colour-matrices corresponds to T x 2 4 = T { . Figure (2) 

shows the diagrams contributing to the colour-ordered partial amplitude Af\q±,g 2 , #3, q~4). 
Squaring the full amplitude A A (qi, g 2 , gs, Q4) one obtains 



A (0) " 



A 
1 



(0) 



(47ra) 2 (47ra s ) 2 |cq| 2 



A^ {qi, g 2 , gz, q 4 ) + Af (q x , g 3 , g 2 , q^ 



|(0), 



o ^ 



1 



(22) 



Note that interference terms between Af'(qi, g 2 , #3, 54) and Af\q\, g 3 , (72, 94) are subleading 
in N c . In four dimensions compact expressions for the partial amplitudes are obtained by 
calculating helicity amplitudes. For the partial amplitude A^' (qi, g 2 , #3, 94) there are only 
three independent helicity configurations, which can be taken to be 



The helicity amplitudes are expressed in term of spinorproducts and read [38]: 



^4 (9i">0a" 5 03":94 ,p 5 >Pe 



-4i 



(4 5) ; 



(12) (2 3) (3 4) (5 6)' 



(23) 



(24) 
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A f ] (if , 9t , 9s ,5 4 >Ps >Pe 



4i 



(31) [12] (4 5) (3|(1 + 2)|6) 



(1 2) S23S123S56 
(34) [4 2] [1 6] <5|(3 + 4)|2> (5|(3 + 4)|2> (3|(1 + 2)|6) 



[3 4] S23S234S56 



(12) [3 4] 523*56 



Af\qt,92,9t^4 ,P 5 >Pet) = 



[1 3] 2 (4 5) (2 1 (1 + 3) |6) + (24) 2 [16] (5|(2 + 4)|3) + [1 3] (2 4) [1 6] (4 5) 



[12] 

S23 s 123 s 56 



(3 4) S23S234S56 



[12] (3 4) 523556 



The notation of spinor products follows [38] . For the singular limits it is sufficient to know 
that spinor products scale as the square root of invariants s^- = ij>i +Pj) 2 '- 



(25) 



For the last two helicity amplitudes there are alternate forms, which have a clearer singular 
structure in the P2IIP3 limit, at the expense of a more obscure behaviour in the triple 
collinear limits pi| j_p2 j \pz and p 2 \ \pz\ \p±- 



A T (Qi, 9t '9 3 , <? 4 , P5 > P& ) = 4 * 



(5 4) [4 2] [12] (5|(3 + 4)|2) 



[2 3] [3 4] Sl23 S234 ( 5 6) 
(3 1) [1 6] (34) (3|(1 + 2)|6> (3|(1 + 2)|6> (5|(3 + 4)|2) 



(26) 



(12) (2 3) Sl23 S234 [5 6] 



(12) [3 4]S 123 S234 



4 0) ( 



11 1 



4i 



g 2 i93i<i4iP^pi) = 

[13] 2 (45)(5|(2 + 4)|3) (24) 2 [1 6] (2|(1 + 3)|6) [1 3] (2 4) [1 6] (4 5) 

[12][2 3]S 12 3S234 (5 6) (2 3)(3 4) Sl 2 3 S 2 34 [5 6] [12](3 4) Sl2 3S234 



As can be seen from these explicit formulae there are simultaneous singularities in single 
collinear limits — ■> and triple collinear limits — > 0. These are nested singularities. 
It should be noted that the helicity amplitudes are sufficient for the derivation of the 
subtraction terms when working in a four-dimensional regularization scheme [39]- [41]. In 
conventional dimensional regularization [42] the squared matrix elements has 0(e)-terms, 
which also have to be taken into account. 



3.2 One-loop amplitudes with one unresolved parton 

For a NNLO analysis of e + e~ — > 2 jets one needs also the amplitude for e + e~ — > qgq up 
to one loop. The colour decomposition of the tree amplitudes and the one-loop amplitudes 
are 

Af ] (91,02,93) = e 2 gc {T 2 ) 13 Af ) (q 1 ,g 2 ,q 3 ), 

A^( qi ,g 2 ,qs) = eVo(T 2 ) 13 (^) 4^(^,02,93). (27) 
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Diagram 2a 




Diagram 2b 



Diagram 2c 



Figure 3: Examples of one- loop diagrams, contributing to different primitive amplitudes: 



Diagram 2a contributes to Ag ^ , 



contributes to A 



diagram 2b contributes to A 



and diagram 2c 



For both the tree amplitude and the one-loop amplitude there is only one partial amplitude. 
However, to obtain a clear factorization in singular kinematical limits, it is convenient to 
decompose the one-loop partial amplitude further into primitive amplitudes [31]. Primitive 
amplitudes can be defined as the sum of all Feynman diagrams with a fixed cyclic ordering of 
the QCD partons and a definite routing of the external fermion lines through the diagram. 
For the case at hand we have 

\NA )M %^) - ^4 )M %U9^) + \N f A^%^ 92Ai ). (28) 

The superscript L or R indicates whether the external fermion line turns left or right 
when entering the loop, while the superscripts [1] and [1/2] indicate the spin of the particle 
circulating in the loop. Examples of one-loop diagrams corresponding to different primitive 
amplitudes are shown in fig. (3). Since we are here only interested in the leading-colour 
contribution, only the primitive amplitude A3 (qi, gi, qs) is relevant. In later sections I 
will drop the superscript "L, [1]" for the leading-colour primitive amplitude. There is only 
one independent helicity configuration, which we can take to be 

Ox '92,%, Pi, Pb- ( 29 ) 
The helicity amplitudes are expressed in terms of spinorpro ducts and read: 

^tf,i},s,K,ri) = 2 ^ (12) ( ( 3 2 4 3) (54) . 

4 >tL ' K (it,a},!h,P4,pi) = vf>4°>(,+ 9 + ?3 -,p 4 -,p+) 

+2 ^<13)[15](34) ^<18)W(45> 

(12)(23) S4 j (12X23)4; v ' 



The renormalized coefficients F , and F\ are given in the HV scheme by 

V 3 (1) = — 1+ (ln(xi) +ln(x 2 ) - ^ J - -R{x 1 ,x 2 ) - ^ln 2 (xi) - ^ln 2 (x 2 ) + \ix 2 
e z \ 6 I e I I 
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3 7 
+- ln(x2) — - + 0(e) + imaginary parts, 

1 (1 - x 2 ) + ln(x 2 ) 

F i = o ^ ^ +0(£J- (31) 

2 (l-x 2 ) 

Here, Xi = S12/S123, x 2 = S23/S123 and the symmetric function R(x\,x 2 ) is defined by 
R(x u x 2 ) = Q ln(xi) ln(x 2 ) - ln(xi) ln(l - x^ + ^2 - ^2(^1)^ + (xi <-> ^2) • (32) 

The amplitude in the FDH scheme is obtained by replacing by 

V 3 % H = V^ + \. (33) 

3.3 The two-loop amplitude 

For a NNLO calculation one also needs the amplitudes for e + e~ — > qq up to two loops. The 
colour decompositions of the tree amplitudes, the one-loop amplitudes and the two-loop 
amplitudes are 

A. 2 \qi,q 2 ) = e 2 c 5 12 A i 2 °\q u q 2 ), 
A { 2 1] (qi,q2) = e 2 c S 12 (^j A^ HM ( qi ,q 2 ), 

A?\ qi ,q 2 ) = e 2 Co^(^) 2 4t M (?i,?2). (34) 

There is only one independent helicity configuration, which we can take to be 

qt,<h, Pa ',Pt- ( 35 ) 
The helicity amplitude expressed in terms of spinorproducts reads: 

4 0) (^? 2 -,P3-,P4 + ) = 2iMM. (36 ) 

S12 

The leading colour contributions to and A 2 2 ^ are: 



A ( 2 l Ua M,Q2,Ps,Pt) = lN c V 2 {1 U¥\qt,q 2 ,Pz,pi)+0(^ 
A?UM^2,P3,Pt) = \(N 2 c -l)vi 2) Af\qtrq2,P3,Pi) + 0(^. (37) 



The leading-colour primitive amplitudes are therefore 

4 1} = V 2 W Af\ 

A? = V^Af. (38) 
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The renormalized coefficients and are given in the HV scheme by [43]- [45] 



v 2 — 




1 


3 




~ 2e 




' 1 

e 


V? = 





7 + f -8 + I* 2 + ksl ^ + ("16 + ^ 2 + ^3 - t^U 4 ) e 2 



12 



3 
2 



-4 + l.) e+ (_ 8+ |,» + ^ 



3 2 

+ 0(e 3 ), 



1440 



433 



J_ 17_ 

2? + 4^3 + ^72"" 

26 . 263 , ^. , 



+ 



4045 83 



1 9083 2153 



432 24^ + 6^ 3 / : 2592 132 
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imaginary parts. 



(39) 



Note that for a NNLO calculation of e + e — > 2 jets one needs the real part of to 
order 0(e 2 ) and the imaginary part of to order 0(e). 



3.4 Ultraviolet renormalization 

The amplitudes in the examples above are the renormalized ones, i.e. the ultraviolet 
subtraction has been performed. To obtain the renormalized amplitudes in the MS scheme 
from the bare ones, one replaces the bare coupling a® with the renormalized coupling a s (/j, 2 ) 
evaluated at the renormalization scale jj 2 : 



1 ,.2e 



a = a s S £ /i 



where 



(4tt) £ e 



(40) 



(41) 



is the typical phase-space volume factor in D = 4 — 2e dimensions, 7^ is Euler's constant, 
and (3q and f3\ are the first two coeffcients of the QCD /9-function: 



Pi = ~ \c A T R N f - C F T R N f , 



with the color factors 



C A = N, C F 



N 2 



1 



■R 



2N ' " 2 

Let the expansion in the strong coupling of the unrenormalized amplitude be 



(42) 



(43) 



11 — z 

AnMre = (4vT«) (47Ta ) 2 C C 



AO) (Oo\ Al 

n,bare \2t{ ) n < 



) 

bare 



+ 



«0\ 2 A2) 



2tt/ 



n,bare 



+ 0(al 



+ subleading colour structures. 



(44) 
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where C is a colour factor. Then, the renormalized two-loop amplitude can be expessed as 
A = 

(47ra) (Aira s )^ c C (S £ ~V 2£ ) 



n — 2 
2 



+ subleading colour structures. (45) 

The relations between the renormalized and the bare amplitudes are given by 

4(0) _ AO) 

n,ren n,barei 



Al) c-l 2e^(l) (w-2) A) ,(0) 

-^n, ren ^£ A* ^n,bare ^ n,barei 

^n.ren ~ J s P ^n.ftare 2 E ^n,bare^ ^ I 4 £ 2 2e n ' bare ' 



(46) 



Below I will discuss the factorization of one-loop amplitudes in singular (soft and collinear) 
limits. The general structure is 

4 1} = SingW-^ + SingW.^l, (47) 

where the function Sing corresponds either to the eikonal factor (soft limit) or a splitting 
function (collinear limit). Eq. (47) holds both for bare and renormalized quantitites. The 
relation between the renormalized and the bare expression for the singular function is 

Smg« = SfV* Sing« e - | SingJl. (48) 

4 Singular limits 

In this section I review the behaviour of amplitudes in singular limits, e.g. when two or 
three partons become degenerate. The knowledge of the singular limits serves as a starting 
point for the construction of the subtraction terms. Most of the material in this section 
is summarized from the literature [19]- [33], although some results, like the triple collinear 
splitting function for colour-ordered amplitudes, are new. When considering singular limits 
it is convenient to work in a physical gauge. In this case, only diagrams where the splitting 
occurs at external lines contribute [46, 24]. In this paper I use the axial gauge n ■ A — 0. 
In this gauge the gluon propagator reads 

—d" u (k,n) = —[- g ^ + n 2 



k 2 k 2 \ k-n (k-n) 

w _ k*rf \ ( v _ V- 
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4.1 Soft gluons 

If a single gluon becomes soft, the partial tree amplitude factorizes according to 

A^( Pl ,p 2 ,p 3 ,...) = Eik^ 1 \p 1 ,p 2 ,p 3 )A^_ 1 (p 1 ,p 3 ,...), 
where the eikonal factor is given by 

S12S23 

and F pfl(T (p) is defined by 

F p ^(p) = g ptM p a -p p g pa . 
The square of the eikonal factor is given by 



Eik(°' 1 )(p 1 ,p 2 ,p 3 ) 



^i3_ 



(50) 



(51) 



(52) 



(53) 



In a similar way the partial tree amplitude factorizes, when two gluons become soft [20] 



with 



4?(Pl,P2,P3,P4,-) = Eik( ' 2 )(p 1 ,p 2 ,p3,p 4 )^ 2 2 (pi,P4,-) ) 

~P\F piia (P2)F aUT (p 3 )pl p{F m(J (p 2 )F aUT (p 3 )p\ 



(54) 



Eik( ' 2 )( Pl ,p 2 ,p 3 ,P4) = 8 



S12S23S34 



P4F p ^(p 2 )F aUT (p 3 )p T 4 

S23 S 34(S24 + s 34) 

The square of the double soft eikonal factor is given by 



S12S23(S12 + Sl3) 

e"(P2)e u ( P3 ). 



(55) 



Eik(°< 2 )( Pl , p 2 ,p 3 ,p 4 ) 
+ — 



(1-pe) 



(S123S34 + S12S234 - s 123 s 2U f 



'14 



+ 



Sl4 



1 1 
+ + 



-2 „2 „2 

3*: 
1 



S 123 S 23 S 234 



Sl2S34Sl23 s 234 S23 \ S 12 S 34 Si 2 S234 Sl23 s 34 S123S234 



(56) 



I introduced the parameter p, which specifies the variant of dimensional regularization: 
p = 1 for the CDR/nV schemes and p = for the FD scheme. In the soft-gluon limit, a 
one-loop primitive amplitude factorizes according to 

A^( Pl ,p 2 ,p 3 ,...) = Eik( ' 1 )(p 1 ,p 2 ,p 3 )-4-i(Pi^3,-) 

+Eik^\ Pl , p 2 ,p 3 ) -A^Upi,P3,-), (57) 
For the example considered here, the relevant one-loop eikonal function is given by [31] 
Eik^ (p 1 ,p 2 ,p 3 ) = 



.^r(i + e )r(i- e ) 



/A~ g 13) 

(-S12X-S23) 
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6e 



Eik^\ Pl , p 2 ,p 3 ). (58) 



Here, 



As in the tree-level case with one soft gluon, the eikonal factor is independent of the variant 
of the dimensional regularization scheme used. 

4.2 Collinear particles 

In the collinear limit tree amplitudes factorize according to 

4 0) (-,P^,...) = ^spiit^^p.o^iC-.P,-)- (so) 

A 

where pi and pj are the momenta of two adjacent legs and the sum is over all polarizations. 
In the collinear limit we parametrize the momenta of the partons % and j as [15] 

. k 2 n 
Pi = zp + k± 



z 2pn' 

Pi = (i _„),_*,-_£_ JL. (61) 

Here n is a massless four-vector and the transverse component k± satisfies 2pk± = 2nk± = 
0. The collinear limits occurs for k\ — > 0. The splitting amplitudes Split*- ' 1 "* are universal, 
they depend only on the two momenta becoming collinear, and not upon the specific 
amplitude under consideration. The splitting functions Split*- ' 1 -* are given by 



Splits = 


7* 


Splits = 




Splits = 





u{p 3 ). (62) 

Due to the sum over spins in eq. (60), spin correlations are retained in the collinear limit 
of squared tree amplitudes. With 

p(o,i) = J2u{p) Split^Split^uft?) for quarks, 

A, A' 

P (o,i) = J2^(P)* Split (0 ' 1) *Split (0 ' 1) £ l/ (p) forgluons, (63) 

A, A' 

and the paramterization eq. (61) one finds 

2 



p(o.i) = —A 



+ (l- p£ )(l-z) 



1-z 
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p(0,l) 



9^99 



3(0,1) 
9^11 



i J 



-5 



/if 



2z 2(1 ^ . „ 

+ J — 4(1 — pe)z(l - z)^j± 



1 — z 



k 



± 



Y v + 4*(1 - ^) 



(64) 



One-loop primitive amplitudes factorize according to [27]-[31] 

A£\..., Pi , Pj ,...) = 

Split(°' 1 )fe,^)4 1 2 1 (...,p,...) + ^ SpUt^(p i ,p i )^i(-,P,-)- (65) 

A A 

Here, the splitting amplitudes Split*- 1 ' 1 -* occur as a new structure. For the one-loop ampli- 
tude A*p for e+e- — > qgq we only need the q — > splitting function. This function is 
given by [30, 31] 



Splits = S^cr 



/i 2 



{h{z) + h) Split£& - / 2 2pi ' 2 e(Pj " ) ti(p i )^ti(p) 



-i Split™ , 
where the coefficients fi(z) and / 2 are 



(66) 



fi(z) = -\ 2 F 1 [l,-e,l-e; 



h 



1-z ' 



1 — pe 



2(l- e )(l-2e)' 
The e-expansion of the hypergeometric function is 



(67) 



(68) 



The hypergeometric function can also be written as 



—2 



1-2 



r(l + e)r(l-e) 



1-2 



+ 1 - z £ 2 Fi (e, e, 1 + e; 1 - z) , 

(69) 



which is useful when studying the soft limit. The interference term of Split*- ' 1 -* with Split*- 1 ' 1 -* 
is given by 



P^ g = J2 U ^ (Splits)* Spht^ u{p) + c.c. = 



-Sl2 



AC*) P^l + h ~i> [1 - pe{l -z)]\- -P<M + c.c. (70) 
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Here "c.c." denotes the complex conjugate. 

For tree amplitudes we also have to consider the case, where three partons become collinear 
simultaneously. In the triple collinear limit we parametrize the momenta as follows [24]: 

Pl = Zip + kxi-—^-, (71) 
where n 2 = 0, 2k±ip = 2k^n = 0, and the Zi and satisfy 

The following formulae are useful: 

2piU 2p i n ( k±i k ±j ' 



™1 i ^"1^ '? 

-/"' 2p ijk n 



( — ~ — J , k 2 u = Zi(s jk - (1 -Zi)s ijk ) . (73) 



Here, pij k = p% +Pj +Pk- For the example of the leading N c contributions to e + e~ — > qggq 
we have to consider only the triple collinear limit q — > qgg. In the triple collinear limit the 
tree amplitude factorizes according to 

A^( Pl ,p 2 ,p 3 ,...) = SplitK s > 1 ,p 2 ,p3)4 l 2 (P,-)- (74) 

A 

Fig. (4) shows the Feynman diagrams contributing to the triple collinear limit q — > qgg. 
We are interested in the square of the splitting function Split qJ?l gg , defined by 

p (o,2) = J2u{p) Split (0 ' 2) *Split (0 ' 2) u(p) for quarks. (75) 

A, A' 

For the splitting q — > qgg I find 

p(0,2) = 
q-*Q99 



-i- A 2(1 - p g ) ( " 3 ~ (1 ~ 2 Xl) " 23)2 + 4(1 - pe)^ + (1 - pe) 2 ^ + 3 - Ape + p 2 e 2 

'123 [ ^2^23 X 2 %1 
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+ 



XiX 2 

1 

Xi 



l + (l-ZzY- P 8zl + 1 + (_1 - Z 23 ) 2 



^4$ 



+- 



^23 

2 2 



^3 

1 + (1 - z 3 ) 2 - pez\ 



Z3Z23 Z3 z 2 3 



Z23 



^23 ^3 



+pe(l - p£)(^23 - £3) - 1 + 3pe] 



1 

+— 

£2 



2 l + (l-^ 3 ) 2 -p^| | 1 + (1 - Z23) 2 - pez 2 23 



z 23 



Z3 



4 , Z3 
Ape — 

^23 ^23 



-3(1 - pe)(z 2 3 - z 3 ) + 2 + Ape}} . 



(76) 



Here x\ = S12/S123, ^2 = S23/S123 and 223 = z 2 + z 3 . Note that eq. (76) gives the 
triple collinear splitting function for colour-ordered amplitudes. This implies a fixed cyclic 
ordering for the two gluons. In contrast, the triple collinear splitting functions given in ref. 
[22] and ref. [24] are the ones for the full matrix element squared, e.g. they include a sum 
over the permutations of the two gluons. 



4.3 Factorizable double unresolved contributions 

For the double unresolved contributions there are also two configurations, which have quite 
a simple factorizable structure in the singular limit. These two cases are: a) the emission 
of two independent pairs of collinear particles and b) the emission of a soft gluon together 
with one pair of collinear particles. Let us first consider the case of the emission of two 
independent pairs of collinear particles. I consider the case p 1 \ \p 2 and p 3 | |p 4 . The amplitude 
factorizes according to 

4 0) (Pi,P2,P3,P4) = £ Split^,^) Split^fo.p,) A® 2 (p,q), (77) 

A, A' 

where p — pi +P2 + 0{k\) and q = P3+P4 + 0(k 2 J. For the leading-colour contributions to 

e + e > qggq only the splittings q — > qg and q — > qg are relevant. Squaring the splitting 

functions in eq. (77) one obtains 

p(0,2) = p(0,l) p(°^) (7Q) 

1 q->qg,q^qg 1 q^qg 1 q-^qg- \' U J 

In the case of the emission of a soft gluon together with a pair of collinear partons it is for 
our purposes sufficient to consider the case of a soft gluon (p 3 ) together with a collinear 
splitting q — > qg of the partons p\ and p 2 - I simply quote the singular factor for the squared 
matrix element [22, 24]: 

p(0,2) _ a S P4 p(0 ,l) (7q) 

Sp3"S34 

Here, p = p\ + p 2 + 0(k\) and the singular factor factorizes into an eikonal factor times a 
splitting function. 
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Figure 5: Splitting topologies for NLO subtraction. At each splitting the emitted particle 
is directed towards the spectator. 



5 Subtraction terms 

In this section I first review the subtraction terms for single unresolved contributions. I 
then give the subtraction terms for double unresolved contributions and for contributions 
with one virtual and one real unresolved parton. 

5.1 Subtraction terms for single unresolved contributions 

The approxiamtion term da^+l is given as a sum over dipoles [15]: 



da 



(o,i) 

71+1 



(80) 



topologies T 



For leading colour contributions the sum reduces to combinations where the emitting parton 
i and the emitted parton j are adjacent, and the spectator k is either adjacent to i or j. It is 
also convenient to denote by k' the other parton adjacent to the pair i and j. Therefore the 
cyclic ordering is either k', i,j, k or k, j, i, k' . The relevant topologies for the example of the 
leading iV c contributions to e + e~ — > qggq are shown in fig. (5). A dipol is constructed from 
amputated amplitudes, evaluated with mapped momenta and dipole splitting functions. 
By an amputated amplitude I mean an amplitude where a polarization vector (or external 
spinor for fermions) has been removed. I denote amputated amplitudes by \A). If the 
amputated parton is a gluon, one has 



If the amputated parton is a quark, one has 

W 0) (-,P,..-)> = 











du(p) 



4 0) (- 



(82) 



The dipole subtraction terms are obtained by sandwiching a dipole splitting function in 
between these amputated amplitudes: 



L (o) , 

l n+l 1 



,Pe 



A 



(o) , 

n+l> 



,Pe 
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Here, c denotes an overall prefactor. For the example e + e — > qggq, it is given by 

c = (47ra) 2 (4ira s ) 2 \c \ 2 x^(N 2 c -l)N c . 



(84) 



The dipole splitting functions read 



Vg^MiPhPiiPk) 

where 



2 1 

Sijk y 
2 1 

$ijk y 

2 1 

s%jk y 



-9 



z{l-y) 

2 



— 9^ - 
T 



\-zi\-y) 
4r(l - yf 



-2 + (l-pe)(l-y)(l-*) 



2 )+(l-p £ ) 4r(1 V) \ ^kl 



k^\ k i 



r = 



ZPk'Pe + ^PePs 



(85) 



(86) 



These splitting functions differ by finite terms from the original choice of Catani and 
Seymour. The new set of splitting functions in eq.(85) is introduced for later convenience 
with the NNLO subtraction terms. Although arbitrary non-singular terms may be added 
to the NLO subtraction terms, the NNLO subtraction terms depend on the actual form of 
the NLO subtraction terms. 

The dipole splitting functions in eq. (85) have also the correct singular behaviour in 
the soft limit. The term 1/(1 — z(\ — y)) equals Sjjfc/( s ij + s jk) and combining this term 
with the one from the dipole, where the roles of emitter and spectator are exchanged, one 
obtains 



Sijk 



+ 



Sijk 



Sijk 



Sij (s«j ~\~ Sjk) Sjk (fijj -|- Sjk) SijSjk 



(87) 



which is the eikonal factor in the soft limit. 

For the q — » qg splitting, the quark polarization matrix ft e is included in the definition 
of the dipole splitting function. For the g — > gg and g — * qq splittings, the full splitting 
function is shared between the dipoles with spectators k and k', where k and k! are the 
partons adjacent to the pair in the colour ordered amplitude A^°\..., k', i,j, k, ...). For 
the g — > gg splitting the dipole splitting function is chosen in such a way that the singularity 
when parton j becomes soft resides in V g ^ gg (pi,Pj;pk), where as the singularity when 
parton % becomes soft resides in V g ^ gg (pj,pi;pk')- The amputated amplitudes are evaluated 
with mapped momenta. The mapping of momenta relates a n + 1 parton configuration to 
a n parton configuration: 



n+l— >ra 



(pi,...,p n+ l) -> (pi,...,p^) . 



(88) 



Such a mapping has to satisfy momentum conservation and the on-mass-shell conditions. 
Furthermore it must have the right behaviour in the singular limits. Several choices for 
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such a mapping exist [15, 37]. One possible choice relates three parton momenta of the 
n + 1 parton configuration to two parton momenta of the n parton configuration, while 
leaving the remaining n — 2 parton momenta unchanged. This mapping is given by 

y 

Pe = Pi+Pj-z. Pk, 

1-2/ 

1 . . 

Ps = Pk- (89) 

1-2/ 

Pi, pj and p k are the momenta of the n + 1 parton configuration, and p e and p s are the 
resulting momenta of the of the n parton configuration. The variables y and z are given 
by 

y = ^, z=^^. (90) 

Sijk $ik t Sj k 

It is easily verified that this mapping satisfies momentum conservation 

Pe+Ps = Pi+Pj+Pk, (91) 

and the on-shell conditions 

Pi = 0, p\ = 0. (92) 
The singular limit corresponds to y — > and in this limit we have 

Pe ~ (Pi + Pj) = 0(y), p s -Pk = 0(y). (93) 
It is also useful to introduce the vector k±, defined by 

y 

y 

With this vector and the value of the longitudinal momentum fraction z, the mapping may 
be inverted to yield 

Pi = zp e + k ± + y (1 - z) p s , 
Pj = (1 - z)p e -k ± + yzp s , 

Pk = (1 - y)p s , (95) 

where y is now given by 

y = n ~u\ V (96) 

Z{1 - Z) (Pe+Ps) 

The singular limit occurs now for k± — > 0. The g — > and g ^ qq splittings also involve 
spin correlation through the spin correlation tensor k^_k v ± . Note that the spin correlation 
tensor k±k± is orthogonal to p e and p s : 

2p e k ± = 2p s k ± = 0. (97) 

Further note that the contraction of p e into an amputated amplitude vanishes due to gauge 
invariance. 

The parameter p specifies the variant of dimensional regularization: p — 1 corresponds to 
the CDR/HV schemes and p = to a four-dimensional scheme. 



k ± = (1 - z) Pi - z Pj - (1 - 2z) Y^Pk- (94) 
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Figure 6: Splitting topologies for NNLO subtraction: Two single splittings from two dif- 
ferent branches. At each splitting the emitted particle is directed towards the spectator. 



5.2 Subtraction terms for double unresolved contributions 

In this subsection I give the subtraction terms for double unresolved contributions. To be 
more precise, these subtraction terms have to approximate 

O n+ 2 da^ 2 - O n+1 o da^l (98) 

in all single and double unresolved limits. Eq. (98) is integrable over single unresolved re- 
gions for observables, which vanish on n-parton configurations, e.g. (n + l)-jet observables. 
However, as discussed in the example of eq. (16), for a generic n-jet observable, eq. (98) 
is in general not integrable over single unresolved regions. Therefore, not only the singu- 
larities corresponding to double unresolved limits, but also the singularities corresponding 
to single unresolved limits have to be subtracted. 

The subtraction term da>n' 2 ^ at NNLO is also written as a sum over topologies: 

da^ = v n 2 \n (") 

topologies T 

These topologies should be thought of as pictorial representations of the pole structure of 
the subtraction terms X>i°' 2 ^ (T) and the way the momenta of the n-parton configurations 
are constructed from the momenta of the (n + 2)-parton configurations. The intuitive 
picture of "iterated dipoles" is only justified for the approximations of da^tl- 

The relevant topologies for the leading N c contributions to e + e~ — > qggq are shown 
in fig. (6) to fig. (7). They fall into two classes: two single splittings from two differ- 
ent branches (topologies and T^) and two sequential splittings from the same branch 
(topologies T c - ). It is sufficient to consider topologies Tjp and T c -Te only. The re- 
maining topologies are related by symmetry to the ones mentioned above. The subtraction 
term is given by 

V^\T) = cJ{T) <^ )(...,p( 2 ), ...,pi 2 ), ...)| ^°' 2 )(T) |A(°)(...,pi 2 ), ...,pi 2 ), ...)> rf^ +2 . 

(100) 
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Figure 7: Splitting topologies for NNLO subtraction: Two sequential splittings from the 
same branch. At each splitting the emitted particle is directed towards the spectator. 



J is a prefactor not singular in any limit, which can be used to absorbe Jacobians into 
the definition of the subtraction terms, thus rendering the analytical integration over the 
unresolved phase space simpler. Here, I simply take J = 1. The splitting function p(°> 2 )(T) 
is sandwiched between double-amputated amplitudes for the topologies and , and 

(2) (2) 

is sandwiched between single-amputated amplitudes for the topologies T c v -T^ . It is 
convenient to write 

v m {T ) = V^(T)-P^(T), (101) 

where is an approxiamation to da^+ 2 an -d ^(o f) * s an approxiamation to da^+l ■ The 

notation with super- and subscripts is as follows: 

VI%{T) (102) 



contributes to daS l,k ^ and approximates a term of da^ l ' ,k '\ (with da"' identified as day 
It should be noted that the subtraction terms are not unique. For example, the complete 

(2) i2) 

triple collinear splitting function P q _^ qgg is shared between the topologies T c , T d and 
Te 2 \ Therefore, a term of the form 

(103) 



2 ' 

S 123 



which is only singular in the triple collinear limit and finite in all other cases, can be freely 
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moved between these topologies. On the other hand, a term of the form 



s 12 s 123 



(104) 



which is singular in the triple collinear limit and in the single unresolved limit S12 — * is 

(2,) 

uniquely associated with topology Tc , as the pictorial representation suggests. Below I 
will give one possible set for the subtraction terms. 

I briefly mention how the n-parton momentum configuration is constructed from the 

( 2 ) 

(n + 2)-parton configuration. As an example I discuss topology T c . For the outermost 
splitting, particle 1 is the emitter, particle 2 is the emitted soft or collinear particle and 
particle 3 is the spectator. Particle 4 is not involved in this splitting. From this 4-parton 
momentum configuration (pi,P2,P3,P4) a 3-parton momentum configuration (pe~\p^\p4) 
is first constructed, using the Catani-Seymour reconstruction function given in eq. (89). 
From the 3-parton momentum configuration a 2-parton momentum configuration (p¥\p^) 
is then constructed, using the symmetric reconstruction functions given by Kosower [37]. 



This symmetric choice is necessary to ensure that topologies and T^' are evaluated 
with the same n-parton momentum configuration. This is necessary for a cancellation of 
(subleading) single unresolved singularities in s± 2 — > and is discussed in detail in sect. 6. 

To present the splitting functions for double unresolved configurations I use the follow- 
ing notation: 



,(2) 



y = i-y, 



l-z, 



(105) 



and pij = pi +pj. 



5.2.1 Topology A 



This topology corresponds to two independent single splittings and the splitting function 
V is inserted into a double-amputated amplitude. The splitting functions read: 



T>(0,2) 

(0,0) q^qg,q^qg 



X 



+ 



1 + 



yizi' 



u 



'13 

2 



^ S 1234 

1 1 



1 



yixi + (2/1 + yizi) x 2 

u 



+ 



+ 



7/1X2 + u 7/1 + yizx 



[-2 + (1 - pe)x 2 ] 



+ 



7/lXi + X 2 7/lXi + 7/1^X2 



-pe (7/1X1 + 7/1X2) 



yi + yizi 



[-2 + (1 - pe)u] + [-2 + (1 - pe)x 2 ] [-2 + (1 - pe)u] 
- 2 + (1 - pe^zi 



(0,2) 

(0,1) q^qg,q^qg 



X- 



1 

y\ 



1-^1(1-7/1) 



ij S 1234 X 1 X 2 L X 1 + X 2 

-2 + (l-pe)(l-y 1 )(l-z 1 ) 



— 2 + (1 - pe) x 2 — pexi 



(106) 
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As already noted, this is inserted into a double-amputated amplitude. The variables xi, 
X2 and x 3 are given by 



xi = y 2 , x 2 = (1 - y 2 )(l - z 2 ), x 3 = (l-y 2 )z 2 . 
The variables y±, zi, y 2 and z 2 are given by 



(107) 



S34 

yi = — , zi 



s 2 4 



y2 



S234 
S12 



'1 - yi) S1234' 

The remaining variable u is given by 



S23 + S24 
z 2 = 1 



•S234 



1 - 1/2) S1234 ' 



(108) 



u 
w 

c 2 



(zi + yiZi) x x + y x z x x 2 + 21X3 + 2(1 - 2w)^/y 1 z 1 z 1 x 1 x 3 , 

i(l-c 2 ), 

S 12 s 23 + S 24 (( s 23 + S24)(Sl3 + S14) — S12S23) — s 14( s 23 + s 24) 



2S23\/ S12S24 (( s 23 + S24)(Sl3 + s 14) — Sl2"S23) 

The amputated amplitude is evaluated with the mapped momenta 



(109) 



l + q + 
1-2- 



y 2 z 2 (1 + q- 2r) 

1 - 2/2^2 

y 2 z 2 (1 + q - 2r) 

1 - 2/2^2 



Pi +rpy + 



1-9 + 



y 2 (l-q- 2r) 



Pi + (1 - r)pW + - 



1 + g 



1-1/2 
_ 2/2 (1 — g — 2r) 
1 - 2/2 



where 



= P34 



1-2/1 

and the variables r and g are given by 

2/2^2 



yi (i) 

P2, Ps ' 



, p2, 
1-2/1 



(110) 



(111) 



2/2 + 2/2^2 

5.2.2 Topology C 

The splitting functions read: 

4 1 



1 + 4r(l - r)y 2 



1-^2 

^2 



;ii2) 



(0,2) 

(0,0) g^ggs 
t 



K 2) ) = ^ 2 



51234 ^1 I 2/1 (2/I+2/124) 



+ 



(j/i + yiz x ) x l + j/!a; 2 j/jxi + u 



, M -- U , —- (-2 + (1 - pe)u) + ^ (-2 + (1 - pefe) 
(2/1+2/1^1)^1+2/1^2 « 
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+- 



y_i 

u 



-77 % 1 + " (- 2 + (! - P £ ) u ) ~ ~ + (! - P £ ) — 

u [{yi + V\Z\) xi + yix 2 \ xi + y x x 2 u u 



(-2 + (1 - pe)y 1 x 2 ) + pe(l - pe)u - pe(l - pe)y 1 x 2 - 1 + 3pe] 



+ (l-p.) 2 ^ + 3-4pe + pV + ^ 
2/i l/i 



-pe 



yi + yi^i 



2 + (1 - ps)y 1 z 1 + - 



-(1 - pe) - 2fc*i + (1 - rfmxi) + u{yi l m) (- 2 + (! " } > 



(0,2) 

(0,1) q-^qgg 

1 



(t c (2) ) = A 



(2). 



4 1 



e o2 



S 1234 ^1 



2 + (1 - pe) x 2 - pexi 



x ■ 



2/i 



2 + (l-pe)(l-j/ 1 )(l-z 1 ) 



l-z 1 (l-y 1 ) 
The variables xi, £2 and x% are given by 

xi=y 2 , x 2 = (1 - y 2 )(l - 22), ^3 = (1 - 2/2)^2- 
The variables yi, zi, 1/2 and ^ 2 are given by 

Sl3 



(113) 
(114) 



yi 
2/2 = 



S12 
S123 

5123 
S1234 



Zl = 



S13 + S23 ' 



Z 2 = 1 - 



1 



S34 



(1 - (l-J/2) S1234' 



(115) 



The remaining variable u is given by 



u 



y 1 z l x 1 + (zi + j/jzj) rr 2 + 21X3 + 2(1 - 2w) v /t7^i^ 2 "x3~, 



c 2 



S 12S23 S 34 + s 13 (( s 13 + S23)( s 14 + s 24) — s 12 s 34) — s 14( s 13 + s 23)' 



2a/si2S 2 3S34Si3 ((Sl3 + S 2 3)(si4 + S24) ~ S12S34) 

The amputated amplitude is evaluated with the mapped momenta 



Pe 2 
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where 



1 



1 + g- 
l-q + 



y 2 z 2 (1 - q - 2r) 

1 - 2/2^2 
y 2 z 2 (1 - q - 2r) 

1 - V2Z2 



rf ) + (l-Opi 1) + - 



l-q 
l + q + 



y 2 (l + q-2r) 



1 - 2/2 

y 2 (l + q- 2r) 
1-2/2 



Pi, 



2/1 (i) 

Pl2 - ~, — P3, Pa 



1-2/1 

and the variables r and q are given by 

2/2 



1-2/1 



P3, 



2/2 + 2/2^2 ' 



g = Wl + 4r(l -r)y 2 



l-z 2 

z 2 



(116) 

P4, 

(117) 
(118) 

(119) 



26 



5.2.3 Topology D 

The splitting functions read: 



(0,2) 

(0,0) q^qgg 
\ 



4 1 



'1234 X l 



yi (yi + yizi) 



(yi + yizi) xi + u 



(yi + yizi) xi + x 2 + yix 3 (2/1 + 2/1*1) xi + u 



- Mi±ig±M! [_ 2 + (1 _ p£) ^ + Ia + Sll3 )] 



+ 



yixi + x 2 + 2/1X3 



[-2 + {\-pe)u\ 



1 

+— 

yi 

-Ape 



+ 



2u 



(1 + + x 2 + I/1X3 xi + u yiXi + x 2 + 2/1X3 



[-2 + (l-peH 



u 



y 1 x 1 + x 2 + 2/1X3 



- (1 - pe) (yixi + x 2 + 2/1X3) + 3(1 - pe)w - 2 + 4pe 



l/i + 3/1*1 



2 + y^i^i 



2(1 -pe) 



y\ [(1 +yi)xi + x 2 + 2/1x3] 
2 



[(1 - 2/1 zif (2/1X1 + x 2 + 2/1X3) - 2m (1 - 2/1Z1) 



+ 



ir 



2/lXi + x 2 + 2/1X3 



^,(0,2) / T (2) 
' (0,1) q-+qgg \ 1 d 



,0) 



6 S 2 



2/1 



l-^i(l-2/i) 



1234 

- 2 



xf 



Xi + x 2 



— 2 + (1 — pe) x 2 — pexi 



+ 8(l-p.)x 3 ! 
xf (xi + x 2 ) 2/1 



+ 



Xi + x 2 



XT 



X2 

2pe— 

Xl 



-(1- 2/1)^1(1-^1) 



The variables Xi, x 2 and x 3 are given by 

xi=y 2 , x 2 = (1 - 2/ 2 )(l - z 2 ), x 3 = (1 - y 2 )z 2 . 
The variables y±, zi, y 2 and z 2 are given by 

•§23 s 13 

2/i = , *i = ■ , 

^123 s 12 + S13 

•Sl23 1 s 14 



2/2 



Sl234 



Z 2 



(1 - 2/l) (1 - 2/2) S1234' 



The remaining variable u is given by 



u = z x x 2 + 2/1*1X3 - 2(1 - 2u?) ^2/1*1*1^2X3, 
w = -(l-c 2 ), 



c 2 = 



s 12 s 23 s 14 + *13 (( s 12 + Sl3)( s 24 + S34) — S^Sm) — Sm( s 12 + S13)' 
2a/ Sl 2 S23 s 13 s 14 (( s 12 + Si3)(s 2 4 + S34) — S^Sm) 
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(120) 



(121) 



(122) 



(123) 



The amputated amplitude is evaluated with the mapped momenta 



Pe 2 



where 



1 



1 + g- 
1 - q + 



y 2 z 2 (1 - q - 2r) 

1 - V2Z2 
y 2 z 2 (1 - g - 2r) 

1 - mz 2 



p(i) + (l_ r )pW + l 



(i) 



l-q 



y 2 (l + q- 2r) 



+ rpW + - 



1 +q + 



1 - 1/2 

y 2 (l + g-2r) 
1-1/2 



P4, 



P4, 



= P23 



l-2/i 

and the variables r and g are given by 

2/2 



1 P1 > 
1-2/1 



2/2 + V2Z2 



1 +4r(l -r)y 2 



1-^2 

Z2 



(124) 



(125) 



(126) 



5.2.4 Topology E 

The splitting functions read: 



?(0,2) 
(0,0) (j^gsg 



(T e (2) ) = # 



1 



x 



S 1234 I ( X l + 2/1^2 + 2/1^3)2/1^2 

4 1 

[-2 + (1 - pe)x 2 ] + 4 - 2(1 - pe)x 2 



+ 



xi + (l + yi) x 2 + yix 3 (yi + y x z x 

' [-2 + (l-pe)x 2 



y\x 2 [(yi + yiZi) x 2 + it] (yi + y ± z ± ) 



Vix 2 



{-pe) 



Vi + Vizi 



-2 + y x z x zx 



(0,1) q-^qgg 

1 



X ■ 



2/1 



\-z x {\-y x ) 



e 2 

S 1234 a 'l a '2 



- 2 



Xi + x 2 



2 + (1 - pe) x 2 - pex 1 



The variables xi, x 2 and x 3 are given by 

a^i=2/2, x 2 = (1 - 2/2) (1 - z 2 ), x 3 = (l-y 2 )z 2 . 
The variables y±, zi, y 2 and z 2 are given by 



(127) 
(128) 



-S23 

2/1 = , zi 



S24 



2/2 



^234 

1 (>Sl 2 3S 2 34 — S23S1234) 



(l-2/l) 



S234Sl 2 34 



S 2 4 + S34 
1 

Z 2 = 1 - 



S234 



(1 -2/2) S1234' 



(129) 
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The remaining variable u is given by 

u = z x x x + yi^xz - 2(1 - 1vo)^/y x z x z x x 1 x z , 
w = -(l-c 2 ), 

S23S34S14 + S24 ((«12 + S13)(S24 + S34) - S23S14) - S12(S24 + S34) 2 on x 

c 2 = , . (130) 

2A/S23S34S14S24 ((Sl2 + Sl3j(S24 + S34) - S23S14) 

The amputated amplitude is evaluated with the mapped momenta 



Pi 2) 



1 + q + 
l-q- 



V2Z2 (1 + q ~ 2r) 

1 - 2/2^2 
y 2 z 2 (l + q- 2r) 

1 - 2/2^2 



Pi + rpy + - 



l-q + 



y 2 (l-q- 2r) 



pi + (1 - r)^ 1} + - 



l + q- 



1-2/2 
3/2 (1 — 9 — 2r) 



1 - 2/2 



131) 



where 



(1) 3/i (1) 
Pe -P23- z Pa, PI ' 



1-2/1 

and the variables r and g are given by 



1-2/1 



P4, 



132) 



2/2^2 



2/2 + 2/2^2 ' 



q = */l + 4r(l -r)y 2 



1 - -22 

^2 



(133) 



5.3 Subtraction terms for contributions with one virtual and one 
real unresolved parton 

In this section I give the subtracttion term dotn ,l \ which approximates 



da, 



+ = (^a) 2 (4na s )^\c \ 2 x](N*-l)N ( 



x 



2vr 

4(0) * , 4 (i) * 4 (o) , 4 (o) * T4 (o) 

^n+l ^n+l T" ^n+l ^Si+l ^n+l "^n+l 



#n+l (134) 



X is obtained from integrating dot£+l over the unresolved phase space and is given for our 
example by 



X = S. 



-1 (4tt) £ |/ 8u 

r(i-e) |U 2 



y(0,l) , v (o,i) 

K <3^93 ^ 9^99 



S23 



Sl2 + S 2 3 



+ 



£23 
/i 2 



S12 + S23 



<3^99 1 9^99 



(135) 
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Figure 8: Splitting topologies for the subtraction terms for the one-loop amplitude with 
one unresolved parton. 



The explicit forms of the functions Vf^qg and Vf-^qg are given in eq. (142). The subtraction 
term dan is written as a sum over topologies: 

dog* = £ V^\T). 



(136) 



topologies T 



The relevant topologies for the leading N c contributions to e + e — > qggq are shown in 
fig. (8). It is sufficient to consider the topology Ta . The topology T^ 1 ' 1 ' is related by 
symmetry to T^^\ The subtraction term is given by 

V^\T) = (137) 

Here, in the first and second term the NLO splitting function p( 0,1 ) is sandwiched between 
An^ and An ■ The third term contains as a new structure the one-loop splitting function 
V^- 1 ' 1 '. It is convenient to write 



V^\T) = V^(T)+V^(T), 



:i38) 



where "P^'q] is (part of ) the approxiamation to da^^ and is an approxiamation to 



da 



(o,i) 

n+l ■ 



(1,0) 

The splitting functions read: 
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Q~*qg ~ g^gg 



y)(i-*) 
*(i-y) 



30 



+ (!-„)- (I-*)" 

The variables y and z are given by 

1/ = 



v (o,i) + y(0,l) 
q^qg ^ g^gg 



Sl2 



l-z(l-y) 



Sl3 



T>(0,1) 



Sl23 



s 13 + s 23 



(139) 



(140) 



The functions g± and #2 are given by 



gi(y,z) = 
1 

~72 



92 



r(l + e)r(l-e) 
1 — pe 



1-z 



1 - (1 - yff 2 F 1 (e, e , 1 + e; (1 - y)(l - z)) 



2(l-e)(l-2e)- 
The sum of the integrated dipoles reads 

Y{-e)Y{-2s)Y{l-ef 



(141) 



v (o,i) + y(o,i) f \ = 
v q^qg 1 y g^ggV I 



r(l - 2e)r(l - 3e) 
4 e(l-e) 



4 + 



8e 



l-3e 



2(l-pe) 



e(l-e) 



3(l-3e)(2-3e) 
The amputated amplitudes are evaluated with the momenta 



(l-3e)(2-3e) 
(142) 



Pe 
Ps 



P1+P2- 



P3- 



1-2/ 



P3, 



1-2/ 



(143) 



For the numerical integration over the (n + l)-parton phase space we need the e-expansion 
of the subtraction terms: 

<4°)|p^)(t)|4 0) > = 4 0) *4 0) 



S12 



+ (3 - 2L) - + 2 Li 2 ((1 -y)(l- z)) + 2 (ln(l - y)f 



+2 ln(l - y) \n(z) + 2 ln(l - y) ln(l - 2) - 2 ln(j/) ln(l - 2) + 2 ln(z) ln(l - 2) 
-2 ln(y) ln(l - y) - \ 2 - ^ ln(l - z) - ^ ln(l - y) - ^ ln(y) + y 



Li-^(i-y) 
+0(e). 



2 + (l-j/)(l-z) 







+ 


— -3 + 2L 







(144) 
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Here, L = ln(si 2 3//i 2 ). The e-expansion of the insertion of V^ 0,1 ^ into the interference term 
with the loop amplitude is given by 

(A^\V^\T) |4 1} > + (AM\V^(T) |4 0) > = 4 0) *4 0) x - 

Sl2 



X 



~ + (-3 + 2L) - - 8 + 3L - L 2 + Jtt 2 
e 2 e 6 



2 + (l-y)(l-z) 



l-z(l-y) 

p(l-y)(l- z )\ + 0{e). (145) 



- + 3-2L 

e 

We see that the poles in e cancel in the sum of the two contributions, as it should be. 

6 Cancellations 

In this section I discuss how singularities cancel between different terms for the various 
singular limits. The most intricate part is the one which is integrated over the n + 2 parton 
phase space: 

J (a+2 dcr% - O n+l o dc£il - O n o da^) , (146) 

where 

dc£ii= E vi n+hn da^= Yl vi n 2 \n ( 147 ) 

rp rp\X) T-i(l) rp(2) 

a )■■■! ^ a )•••! ^ 

Although individual terms in the integrand of eq. (146) are singular, the singularities cancel 
in the sum and the integrand is integrable. Therefore the integration can be performed by 
Monte Carlo methods. It is instructive to see in detail, how the cancellations of singularities 
take place. Critical region in phase space are regions where one or two particles are 
unresolved. Double unresolved contributions can be divided into the following classes: 

a) Two pairs of separately collinear particles. 

b) Three particles collinear. 

c) Two particles collinear and a third soft particle. 

d) Two soft particles. 

e) Coplanar degeneracy. 

The last item occurs only in the subtraction terms and is an artefact of our choice for 
the NLO subtraction terms. It is included here for completeness. The single unresolved 
contributions are divided into two classes: 
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f) Two particles collinear. 

g) One soft particle. 

For the leading-colour contributions to e + e~ — > qggq it is sufficient to check the following 
cases: Particles f and 2 collinear and particles 3 and 4 collinear (case a); particles f, 2 
and 3 collinear (case b); particles f and 2 collinear and particle 3 soft (case c); particles 
2 and 3 soft (case d); the sum of the four-momenta of particles 2 and 3 linear dependent 
with the four-momenta of particles 1 and 4 (case e); particles 1 and 2 collinear (case f); 
particles 2 and 3 collinear (also case f); particle 2 soft (case g). All other cases are related 
by symmetry to the ones above. 

All limits have been checked using the symbolic algebraic manipulation program "FORM" 
[47]. 

6.1 Double unresolved regions 

In the double unresolved cases, the singularities in the NLO subtraction terms da^+i are 
cancelled by the subtraction terms dct^ 2 } . In particular we always have 



O n +l 






7,(0,2) 
(0,1) 


n(?f) 


-o n 


7,(0,2) 
(0,1) 




= integrable, 




O n +l 


-7,(0,1) /^(lK 
u n+l\ 1 b ) 


-o n 


7,(0,2) 
(0,1) 




-o n 


7,(0,2) 
(0,1) 


„(?f) 


= integrable, 




On+1 




-o n 


7,(0,2) 
^(0,1) 




-o n 


7,(0,2) 




= integrable, 




On+1 


Wo.iWiK 

u n+l\ 1 d > 


- o n 


7,(0,2) 
^(0,1) 




- On 


7,(0,2) 
^(0,1) 


ST?) 


= integrable, 


(148) 



for all double unresolved cases. Here, "integrable" stands for "integrable over the double 
unresolved region" . It remains to discuss the cancellations between da n °l 2 and da^'^ n in 
the double unresolved regions. It should be noted that in specific double unresolved limits 
some subtraction terms in eq. (148) may be integrable by themselves. This will also be 
briefly discussed below case by case. 

6.1.1 Case 1 and 2 collinear and 3 and 4 collinear. 

This case corresponds to the emission of two independent pairs of collinear particles (pi | \p 2 
and P3IIJ94). For each pair of collinear particles we write the corresponding momenta 
according to eq. (61). We are interested in contributions, which scale in the double 
collinear limit as |/cj_| -4 . Cancellations of these terms occurs in the following combination: 

O n+2 da n % - On n (Ti 2 )) - O n V$f } n (Tf ) = integrable. (149) 
For completeness I also state the cancellations for terms related to da n °^ : 
O n+1 V^(T^)-O n V^f )n (T^) = integrable, 
O n+1 V^(T^)-O n V^ )n (T^) = integrable. (150) 
All other subtraction terms are separately integrable in this limit. 
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6.1.2 Case 1, 2 and 3 collinear. 

This case corresponds to the triple collinear limit (p±\ \p 2 \ \pz)- We parameterize the mo- 
menta of the three collinear particles according to eq. (71). We are interested in contri- 
butions, which scale in the double collinear limit as |/c±| -4 . Cancellations of these terms 
occurs in the following combination: 

O n+2 da<® 2 - O n n (Tf) - O n B (Tf ) - O n jTf) = integrable. 

(151) 

For completeness I also state the cancellations for terms related to daf^ '• 



l ra+l 



O n+ iV^(TU)-O n V§f )n (TP) = integrable 



On + i^(T^)-O n v\° f,(TP) = integrable, 



(0,1) n\ 

O n+1 V^(T^)-O n V^f )n (T^) = integrable. (152) 
All other subtraction terms are separately integrable in this limit. 

6.1.3 Case 1 and 2 collinear and 3 soft. 

This case corresponds to the emission of a collinear pair (p 1 1 \p 2 ) together with a soft gluon 
(p 3 ). In this limit, the invariants si 2 , s 23 , S34 and S123 become small. We are interested in 
contributions, which have three of these four invariants in the denominator. Cancellations 
of these terms occurs in the following combination: 

O n+2 daf +2 - ^( ( o°o)n( T ) = integrable. (153) 



71 T"i(2) 

J J a T"i ± g 



For completeness I also state the cancellations for terms related to dot£+l'- 

O n+ i V^l{TV) - O n P^JjTf ) - O n V^] n (T^) = integrable, 
O n+1 V^(TW) - O n V§f )n (T^) = integrable, 

)(0,2) 
(0,1) 

All other subtraction terms are separately integrable in this limit. 



On+i ^fiiW) - O n <'?J n (Ti 2) ) = integrable. (154) 



6.1.4 Case 2 and 3 soft. 

This case corresponds to the emission of two soft gluons (with momenta p 2 and p 3 ). If we 
scale the momenta of two soft particles as 

p 2 -> \p 2 , p 3 -> Ap 3 , (155) 



34 



we are interested in contributions, which scale in the double soft limit as A 4 . Cancellations 
of these terms occurs in the following combination: 

O n+2 da%- °n V tmn( T ) = integrable. (156) 

c >—i 1 h 

For completeness I also state the cancellations for terms related to dot£+l ■ 

On+i VSS(TP) - O n V%% n {Tf ] ) - O n V§f } n (lf ) = integrable, 

,(0,2) rT (2)x _ D v (0,2) W 
(0,1) n y 1 e ) ^(0,1) n y 1 9 

All other subtraction terms are separately integrable in this limit. 



O n+ i V^}(TV) - O n rfg „(Ti 2 )) - O n V%% B (T«) = integrable. (157) 



6.1.5 Coplanar case 

For completeness I briefly discuss singularities occuring in coplanar configurations. The 
matrix element squared da^_ 2 ^ s n °t singular in this limit. This singularity occurs only 
in the subtraction terms. An example is given by the subtraction term V^^T^). This 
subtraction term has a singularity in 

2p 4 p« = 51235234 = 52351234 . (158) 

Sl23 — s 23 

The numerator is a Gram determinant and can also be written as 

S123S234 - S23S1234 = ( 1 - 1 2 + 3 14- ) (4 - 1 2 + 3 1 1 - ) . (159) 

This expression vanishes, whenever the sum of the four-momenta p 2 + P3 becomes linear 
dependent with pi and p 4 . Therefore, 2p^pi > becomes small, whenever this condition is 
fullfilled and whenever p 2 + p% is not linear dependent with p 1 . (In the later case also the 
denominator in eq. (158) vanishes.) This type of singularity occurs only in the subtraction 
terms and the cancellation takes place in the following combination: 

O n+1 V^}{T^)-O n V^ n {T^) = integrable. (160) 



6.2 Single unresolved regions 

The expression in eq. (148) is integrated over the complete (n + 2)-parton phase space. 
This implies in particular an integration over regions which correspond to single unresolved 
configurations. It is clear that in this case we cannot expect a cancellation of singularities 
between terms, which involve O n+2 and O n , e.g. similar to that what we had in the double 
unresolved case. In the single unresolved case the (n + 2)-parton configuration used in 
the calculation of the observable O n+2 will go smoothly to a (n + l)-parton configuration. 
The numerical value for the observable O evaluated with this (n + l)-parton configuration 
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will however in general have no relationship with the numerical value obtained from the 
evaluation with the n-parton configuration of the subtraction terms. A similar argument 
applies to the factorization of the matrix element. Both effects spoil a cancellation of the 
singularities. The cancellation of the singularities has to occur therefore either between 
terms involving O n+ 2 and O n +i, or between terms involving O n+ \ and O n , or between 
terms involving only O n . In the following, "integrable" stands for "integrable over the 
single unresolved region". 



6.2.1 Case 1 and 2 collinear. 

This case corresponds to the collinear limit Pi||p2- We parameterize the momenta of the 
two collinear particles according to eq. (61). We are interested in contributions, which 
scale in the collinear limit as |A:_l| -2 . A cancellations of these terms occurs in the following 
combination: 

O n+2 da^ 2 - O n+1 V™(tP) = integrable. (161) 

This is just the combination which already occurs at NLO. At NNLO we also have to take 
into account the singularity in V^+l(T^). This one is cancelled in the combination 

O n+1 l?£.V(rW) - O n V§f )n (T?>) = integrable. (162) 

Finally, there is a third equation, describing the cancellations among the NNLO subtraction 
terms: 

On 2>gg n (Tf ) + O n V^ n (T? ] ) ~ O n PgjjTf ) - O n V^] n (TP) = integrable. 

(163) 

All other subtraction terms are separately integrable in this limit. Note that the pieces 

On < 2 j „(Tf ) - O n V^] n (Tf) (164) 

and 

On 2>g M 2) ) - On V^] n (Ti 2 )) (165) 

are with our choice of subtraction terms in general not integrable separately. Only the 
sum of the two pieces is integrable. In order for the cancellations to occur, the n-parton 
configuration of topology must approach in the limit S\2 — > the same configuration as 
the n-parton configuration of topology T c . With a symmetric choice for the reconstruction 
functions in eq. (110) and eq. (117) the n-parton configurations for topologies and 
are identical for any value of s 12 and this requirement is fullfilled trivially. 
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6.2.2 Case 2 and 3 collinear. 



Here we discuss the case, where particles 2 and 3 become collinear. Again, we parameterize 
the momenta of the two collinear particles according to eq. (61). We are interested 
in contributions, which scale in the collinear limit as |/c±| -2 . The NLO relation for the 
cancellation of singularities reads: 

O n+2 da% - O n+1 V^(T^) - O n+1 vS#(Ti») = integrable. (166) 

In addition, there is one relation among the subtraction terms with n-parton configurations: 

d t^ ' 2 ) /t^( 2 )\ I /n f)(0> 2 ) ( r r(' 2 ) \ _ /n t-)(0> 2 ) ( r r(. 2 )\ _ (n T)( > 2 ) /' r r( 2 ) N \ 

Wn U {fifi)n\ 1 d ) ^~ w n ^(0,0) n^ 1 h ) w n ^(0,1) n \ 1 d ) W n ^(0,1) n\ L h ) 

+O n B (Tj 2 >) + O n V™ B (T< 2 >) - O n V%* „(Ti 2 )) - O n V™ n (lf) 
= integrable. (167) 

All other subtraction terms are separately integrable in this limit. Again we have to require 
that the n-parton configurations of the topologies Te 2 \ Tg 2 ^ and all approach the 
same n-parton configuration in the limit S23 — > 0. With the choice for the reconstruction 
functions in eq. (124) and eq. (131) this requirement is fullfilled. 

6.2.3 Case 2 soft. 

This case corresponds to the emission of a single soft gluon (particle 2). If we scale the 
momentum of the soft gluon as 

P2 -> Ap 2 , (168) 

we are interested in contributions, which scale in the soft limit as A~ 2 . Cancellations of 
these terms occurs in the following combination: 

O n+2 da^ 2 - a+i V^liT^) - a+i iffiCzf) = integrable. (169) 

This is the usual NLO relation. In addition, we have 

On V%f )n (Ti*>) + O n PgjTi 2 )) - O n V^ )n (T^) - O n P[ °3 JT] 2 )) = integrable, 

On V%* n (Tf>) + O n PgjTf ) - O n vf^ n (T [ d 2) ) - O n PgjTf ) = integrable. 

(170) 

All other subtraction terms are separately integrable in this limit. 

6.3 Cancellations for contributions with one virtual and one real 
unresolved parton 

The cancellation mechanism for contributions with one virtual and one real unresolved 
parton are quite similar to those already encountered at NLO. The only difference is the 
functional form of the subtraction terms (as given in eq. (139)), which approach smoothly 
the appropriate singular limit (as given in eq. (70)). 
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7 Outlook and conclusions 



In this paper I considered the extension of the subtraction method to next-to-next-to- 
leading order. Such a method is needed to extend fully differential perturbative calcula- 
tions in quantum field theories from NLO to NNLO. As a particular example I considered 
the leading-colour contributions to e + e~ — > 2 jets and derived the subtraction terms at 
NNLO. The subtracted matrix elements can be integrated numerically over the appropriate 
phase space. The method presented here is general and not restricted to the example of 
e + e~ — > 2 jets. Subtraction terms for other splittings (like g — > ggg) and other kinemati- 
cal configuartions (e.g. with partons in the initial state) can be worked out along the same 
lines. 

The subtraction terms still have to be integrated analytically over the unresolved phase 
space. With the advance of integration techniques witnessed in the past years, this seems 
feasible. The completion of this program will open the door to fully differential NNLO 
Monte Carlo programs. 
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